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Abstract. We provide an example of how the complex dynamics of a recently introduced model can be 
understood via a detailed analysis of its associated Riemann surface. Thanks to this geometric description 
an explicit formula for the period of the orbits can be derived, which is shown to depend on the initial 
data and the continued fraction expansion of a simple ratio of the coupling constants of the problem. For 
rational values of this ratio and generic values of the initial data, all orbits are periodic and the system is 
isochronous. For irrational values of the ratio, there exist periodic and quasi-periodic orbits for different 
£N) initial data. Moreover, the dependence of the period on the initial data shows a rich behavior and initial 

, | , data can always be found such the period is arbitrarily high. 
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1. Introduction 



Although the evolution of physical systems takes place over real time and can be described by real variables, 
the interest of extending the study to complex values has long been noticed. For instance, the solutions of 
a nonlinear system of ODEs representing a physical phenomenon might have singularities off the real time 
axis whose position depends on the initial data. Even though the time variable is followed along the real 
axis, the position of these complex singularities provide information on the evolution of the system and its 
dependence on the choice of initial data. 

Historically, one of the first successes of this paradigm is the collection of techniques now known as Painleve 
analysis, originally introduced by Painleve and Kowalevskaya [I] [2]. In essence, they consider an ansatz of the 
local behaviour of a solution near a singularity in terms of a Laurent series, introducing it in the equations 
and determining the leading orders and resonances (terms in the expansion at which arbitrary constants 
appear). Painleve analysis has been extended to test for the presence of algebraic branching (weak Painleve 
property [5]). These analytic techniques (which have been algorithmized and are now available in computer 
| packages [4]) constitute a useful tool in the investigation of integrability: in many nonlinear systems where 

no solution in closed form is known, Painleve analysis provides information on the type of branching featured 
by the general solution or by special classes of solutions. It has also proved useful to identify special values of 
the parameters for which generally chaotic systems such as Henon-Heiles or Lorenz are integrable [6] . The 
idea that obstruction to integrability is encoded in the branching properties of solutions and, more precisely, 
that dense branching is responsible for non integrability, preventing from extracting any useful information 
on existing first integrals, was first put forward by Kruskal in connection with his poly-Painleve test [7J|8], a 
generalization of the Painleve test allowing one to investigate, in particular, the degree of multivaluedness of 
the solutions. 

Tabor and his collaborators initiated the study of chaotic systems from the point of view of the singularity 
structure od their solutions [HI HI UHl HU [H] ■ Their local analytic approach (^/-series) was complemented by 
numerical techniques developed for finding the location of the singularities in complex time and determining 
the order of branching [13] . Similar studies relating singularity structure, chaos and integrability have been 
performed by Bountis and his collaborators. They recognised that local analysis alone would not be able 
to characterise global asymptotic properties of the systems. In a number of papers |14| 115) they propose to 
call integrable those cases in which the Riemann surface has a finite number of sheets, and non-integrable 
if the number is infinite. Using mostly numerical evidence they conjecture that in the non-integrable cases 
the Riemann surfaces are infinitely-sheeted and the projection on the complex plane of the singularities is 
dense. Combining analytical and numerical results for a simple ODE, Bountis and Fokas [16 have identified 
chaotic systems with the property that the singularities of their solutions are dense. More recently, a research 
program to investigate classical dynamics extended to the complex domain is being carried by Bender and 
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his collaborators, [T71 [TH1 UP , 20 . In these works extensive use is made of numerical explorations to describe 
complex and interesting behaviours, which include chaotic systems |21) . but the necessity to study the 
Riemann surfaces is recognised. 

A precise correspondence between dynamical properties of the system (e.g. chaotic behaviour or sensitive 
dependence) and specific geometrical properties of the Riemann surface is still an open problem, and this 
paper can be seen as one more step in this direction. Our motivation was to introduce a model which is simple 
enough that a full description of its Riemann surface can be performed, yet complicated enough to feature a 
rich complex behaviour, whose description would not be feasible by standard techniques of qualitative theory. 
Such a model was initially presented in [22], and in a subsequent paper [23] some results were announced 
without proof since they required a full description of the Riemann surface, a task which is carried out in 
this paper. We stress that achieving the complete description of the Riemann surface and its consequence for 
the determination of the period of the solutions involves the use of various mathematical techniques ranging 
from complex analysis to combinatorics and showing remarkable connections with number theory (continued 
fraction expansions). 

An important element of the construction of the model is a change of dependent and independent vari- 
ables which becomes useful to identify isochronous systems [231 [55] . Using local analysis and numerical 
integration in two many-body systems in the plane [371 12E US] j it was discovered that outside the isochrony 
region there exist periodic solutions with much higher periods [3D], as well as possibly aperiodic solutions, 
and the connection among this phenomenology and the analytic structure of the corresponding solutions as 
functions of complex time was illuminated. The model studied in this paper was introduced as a prototype to 
understand the complex behaviour with many periodic orbits observed numerically in the models described 
above. 

Other attempts to study models which are able to produce chaotic motion and yet lend themselves to a 
complete description of their Riemann surfaces are those related to the inversion of hyperelliptic integrals 
[3T][32]. 

This paper is organized as follows: in Section 2 we present the model, show that it can be reduced to 
quadratures and how the general solution can be written in terms of a multi-valued function that motivates 
the study of an associated Riemann surface. In Section 3 we describe the geometric strucure of the Riemann 
surface, originally for the algebraic case (/i £ Q) . The connection with graph theory and Ferrer diagrams 
is illustrated in Section 4, which shows also how some results can be framed within the theory of continued 
fractions. The formulas that allow the explicit determination of the period are given in Section 5, both for 
the rational and irrational cases. Finally, some conclusions and further work are outlined in Section 6. 

2. The model 

The model we analyze in this paper is given by the following set of three coupled first order ODEs: 

/ -I \ ... 9n+2 g n +l 

(1) z n +iuiz n = — — 1 • 

Notation: Here and hereafter indices such as n, m range from 1 to 3 and are defined MOD (3), where mod 
is defined for all positive integers a and b as follows 



(2) a mod (b) = 



a MOD (b) if a mod (b) ^ 
b if a mod (6) = 



(hence, for instance, 3 mod (3) = 3, 4 mod (3) = 1, the usefulness of this notation will become apparent 
below). The dependent variables z n = z n (t) are complex functions and indicate the positions of three 
interacting bodies in the plane; the independent variable t ("physical time") is real and the superimposed 
dots denote differentiation with respect to t. The parameter ui is strictly positive and is associated to the 
period 

(3) T=l 

We fix hereafter u> — tt without loss of generality so that the fundamental period T = 1. 

The three quantities g n are arbitrary "coupling constants" (possibly also complex; but in this paper we 

restrict consideration only to the case with real coupling constants). 
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In the following we will focus only on the "semisymmetrical case" characterized by the equality of two of 
the three coupling constants, say 

(4) 9i = 92 = 9 , 93 = f, 

since in this case the treatment is simpler yet still adequate to exhibit most aspects of the phenomenology 
we are interested in (see |23j). In this semisymmetrical case it is convenient to introduce the constant /j,, 

(5) M = 7T^ 

whose value, as we shall see, plays an important role in determining the dynamical evolution of our model. 

In [53] the general solution of system ((TJ) was given. In the semisymmetrical case (g]), if \i ^ and /j ^ 1, 
the solution of (JTJ) reads as follows. 

M rjtt 2z 3 (0)-z 1 (0)-z 2 (0) r , 9 . , Ull i/2 
z s (t)=Ze — — [rj exp(~2^7^^) + 1J 

(6a) • (-Mi)] 1/2 + (-iy [12 M -3u,(i)] 1/2 ), a = l,2, 

/fiiN /.\ ry fcrt 2 Z 3 (0) ~ Z\ (0) ~ Z 2 (0) . , . . .,1/2 1/2 

(6b) z 3 (t) = Ze l " T — — [77 exp(-2«7ri) + 1] [w (t)\ 

where 

(7) Z = Zl + Z2 + Zz 

3 

is the center of mass, w (t) is the solution of the nondifferential equation 

(8) [w(t) - Htjf = i? exp(2i7rf) + £ = i? [exp (2 i tt 4) + rj\ , 
and the constants R, £ and 77 are defined in terms of the initial data: 

1 tt { [zi (0) - z 2 (0)] 2 + [z 2 (0) - z 3 (0)] 2 + [z 3 (0) - zi (0)] 2 } 

(9a) " = — HfT^ L - 1 - 

(9b) i?= 3(/ + 8g) ffll-K]"- 1 , 

2*7T [2z 3 (0)-zi(0)-z 2 (0)] 2 

(9c) t = Rv, 

where n in f|9b[) is given by 

(Qd) r _ 2 M [2z 3 (0)-z 1 (0)-z 2 (0)] 2 



[«i(0) - *a(0)] 2 + [z 2 (0) - 2 3 (0)] 2 + [z 3 (0) - zi(0)] ' 



At this stage, it is mandatory to note that f|9b[) contains a degeneracy of order fjt — 1. In the following sections 
we will explain the role of such a degeneracy in the construction of the solution and how to remove it. 
If we now set 

(10a) w{t) = w [f (i)] , 

(10b) £ = R cxp(2i7rt) +£ = i? [exp (2 i 7r t) + 77] , 

we can rephrase © as follows 

(11) Ko-ir -1 Mor"=£- 

Note that this equation is independent of the initial data; it only features the constant /1, which only depends 
on the coupling constants, see (J5J) . Moreover, it defines the Riemann surface T consisting of points (£, w) € T 
such that ([TT]) is satisfied. 

As £ travels in the complex £-plane on the circle E defined by (jlObp , the dependent variable w (£) travels 
on the Riemann surface determined by its dependence on the complex variable £, as entailed by the equation 
(fTTj) that relates w (£) to its argument £ - starting at t = from £ = £0, 

(12a) £o=Z + R=(r 1 + l) R, 
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(see (j9a) - (f9d]) ) and correspondingly from iu(£o) =Wo, 

,, 9W i MO) - * 2 (Q)] 2 + MQ) - z 3 (p)] 2 + [z 3 (Q) - *i(o)] : 

U^DJ Wq — — = — — — — ~2 



k 2/x [2z a (0) -zi (0)-%(0)]' 

We observe that £ and £o feature the same degeneracy as R. 

We still must fix the degeneracy of the square roots appearing in (|6a |) - (|6b|) . Since there is no degeneracy 
in the determination of wo and rj, the determination of the signs of the square roots in (p3a]) - (|6bp is fixed by 
demanding that these formulae are consistent with the initial data at t — 0. 

Despite the periodicity in time of £, the corresponding time evolution of w via (jllj) could be much more 
complicated (possibly aperiodic) due to the multivalucdncss of w as a function of £. This evolution must be 
studied by lifting a circular path to the Riemann surface defined by (JTTJ) . This motivates the necessity to 
study the geometric structure of the Riemann surface, which we address in the following section. 

3. The Riemann surface 

In this section we discuss the structure of the Riemann surface V consisting of points (£, aijcT such that 
(JTTJ) is satisfied. 

For rational values of y, = p/q, which is the case if both the coupling constants / and g are rational numbers, 
see (O, equation describes an algebraic curve. We treat the case of irrational /i via an appropriate limit 
of the case with rational fi. Moreover, as already pointed out in [23] and as we recall in the next section, one 
needs to distinguish among the three cases with fi < 0, < fj, < 1 and fi > 1. 

We start the analysis of the structure of T by ignoring the dependence of £ on the "physical" time t (see 
(jlObp ). and by considering instead £ as an independent complex variable, evolving on a generic (possibly 
closed) path on the complex £-plane; only after having thereby obtained an appropriate understanding of the 
topological properties of T, we proceed and analyze the consequences of the "physical movement" of £ along 
the circle 2, as described by (|10b|> . 

3.1. Movable and fixed singularities. The function w(£) defined by (fTTj) features two types of singulari- 
ties: the "fixed" ones occurring at values of the independent variable £ - and correspondingly of the dependent 
variable w - that can be read directly from the structure of ljll|): and the "movable" ones occurring at values 
of the independent and dependent variables, £ and w, that cannot be directly read from the structure of (jTTJ) 
(they "move" as the initial data are modified). 

In order to investigate the nature of the movable singularities, it is convenient to differentiate (JTTJ) , thereby 
obtaining (by repeated use of (fTi"j) ) 

(13) £w = , 

w — fl 

where the prime indicates differentiation with respect to £. The position of the movable singularities, £b, and 
the corresponding values of the dependent variable, Wf, = w(£b), are then characterized by the vanishing of 
the denominator in the right-hand side of this formula, yielding the relation 

(14a) w b = n, 

which, combined with (fTTj) (at £ = £b) is easily seen to yield 

(14b) a-^-rexp^^fc) , k = 1,2,3,... 

' fi — 1 



(14c) r = (fi - 1) 

V M 

In (|14c[) it is understood that the principal determination is to be taken of the /z-th power appearing in 
the right-hand side. Formula (|14bj) shows clearly that the number of these branch points is infinite if the 
parameter fi is irrational (fj, ^ Q), and that they then sit densely on the circle B in the complex £-plane 
centered at the origin and having radius r, see (|14c[) . On the contrary, if /i is rational (/i G Q) the branch 
points sit again on the circle B in the complex £-plane, but there are only a finite number of them. As proved 
in [23], these movable singularities are all square- root branch points. 

Then, let us consider the "fixed" singularities, which clearly can only occur at £ = oo and at £ = with 
corresponding values for w. 
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Two behaviors of w(£) are possible for £ oo, depending on the value of (the real part of) fj,. The first is 
characterized by the ansatz 

(15a) w(0=af + o(\af) , < 0, 

and its insertion in (jllj) yields 

(15b) /3 = — — , — - exp(i7r/i), 

which is consistent with (|15a[) iff 

(15c) /i>0. 
The second is characterized by the ansatz 

(16a) ^(0 = l + ae 3 +o(|C| /3 ) , /3<0, 

and its insertion in (TlTj) yields 

(16b) /3 = — — r j 

A 4 - 1 

which is consistent with (|16a|) iff 

(16c) ^ < 1 . 

We therefore conclude that there are only three possibilities: 

• if n > 1, only the first ansatz, ()15a,|) - (|15c|) . is applicable, and it characterizes the nature of the branch 
point of w(£) at £ = oo; 

• if fi < 0, only the second ansatz, (|16ap - (116cp . is applicable, and it characterizes the nature of the 
branch point of w(£) at £ = oo; 

• if < fi < 1, both ansatze, (|15ap - (|15cp and (|16ap - (|16cp . are applicable, so both types of branch points 
occur at £ = oo. 

The special cases /i = and /i = 1 require a separate treatment |33j . Next, let us investigate the nature of 
the singularity at £ = 0. Two behaviors are possible: either 

(17a) w(0 = a^ +o{\^f) , > 0, 

(17b) /3 = --, a M = -exp(zTT^) , 

A* 

which is applicable if and only if 

(17c) fj, < ; 

or 

(18a) w (0 =l + a^ + o(|£|") , > 0, 

(18b) /3 = — L- , o^-^l, 

M 1 

which is applicable if and only if 

(18c) a* > 1 . 

This analysis shows that the function u>(£) features a branch point at £ = the nature of which is characterized 
by the relevant exponent (3, see (11 7b[) or (I18bp . whichever is applicable (see (|17c|) and (I18c|) ). There is no 
branch point at all at £ = if neither one of the two inequalities (|17cp and (|18c|) holds, namely if < /i < 1. 

Moreover, we observe that the case /i < can be immediately worked out from the case fi > 1 via the 
following replacement 

(19) wkw-I, jui-^l — /x. 

For this reason, only the two cases < /i < 1 and /i > 1 lead to different Riemann surfaces. Although the 
geometrical properties of the Riemann surface are quite different in the two cases, the techniques employed 
in their study are essentially the same. For this reason, to avoid unnecessary repetitions we concentrate in 
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this paper only on the first case < \x < 1, which, as we shall see, leads to a rich and complex behavior. 
The second case \i > 1 which also has interesting consequences for the dynamics, specially in the case fi ^ Q, 
shall be discussed in a subsequent paper. 

3.2. The case < \i < 1 and fi € Q. If /.i is a rational number in the open interval (0, 1), namely if 

(20) n e Q and < fi = - < 1 , 

q 

where p G N and g 6 N + arc coprime natural numbers, then ([lip implies that T is an algebraic Riemann 
surface characterized by the polynomial equation 

(21) (w~l) q - p w p £ q = 1 with 0<p<q, 

which defines the q-valued function w — ttf(£). Since, V£ G C, the polynomial (|21l) admits q complex roots 
and each root corresponds to a sheet of the Riemann surface T, it follows that T is a g-sheeted covering of 
the complex £-plane. 

Condition (|2"0j). via formulae (|15aj) - (jl5cj) and (|16ap - (|16c|) . entails that, for £ w oo, the oo- configuration of 
the q roots of (|2ip consists of p roots lying on a small circle of radius around the origin, and of 

(q — p) roots lying on a small circle of radius 0(|£| - i-r ) around 1 (see Figure [T]). From the point of view of 
the Riemann surface, we see that, at £ = oo, the branch point (oo, 0), of order (p — 1), connects p sheets and 
the branch point (oo, 1), of order (g — p — 1), connects (q — p) sheets. 



Figure 1. The oo- configuration of the q roots of (|21|) . for q = 12 and p = 5. The labeling 
of the roots is explained in the text. 



In the finite part of T, condition (|20|) and formulae (|14a[) - (|14cj) imply that there are only q square-root 
branch points: 

(22a) K b W ,|i)6r, j = l,...,q, 

defined by the equation: 

(22b) - — 



pP{q-p)i-P 

These q square-root branch points correspond to the collision of a pair of roots of (j2Tj) ; they are clearly 
located on the circle B in the complex £-plane centered at the origin and having radius rb 

(22c) r b = ^(i^f>0. 

q-p V P J 

It is convenient to sort these q square-root branch points sequentially in counterclockwise order (see Figure 
[5]); an opportune choice of the first square- root branch point clearly arbitrary at this stage, is suggested 
by the direct problem and will be discussed later. 

The genus of T is 0; this is an immediate consequence of the Hurwitz formula, V = 2(J + G — 1), where 
V is the ramification index of the surface, J is the number of sheets and G is its genus. In our case: J = q 
and V = q + (q -p - 1) + (p - 1) = 2(q - 1), entailing G = 0. 

Moreover, equation (|2ip exhibits several symmetries. The ones used below are: 
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Figure 2. The cut £-plane of the Riemann surface Y with q = 12 and p = 5 is the interior 
of the oriented contour 7. 

• Symmetry under a 2?t/q rotation. The set of roots is left invariant by a rotation around the origin of 
the £-plane by an angle 27r/g (and of course any integer multiple of it). 

• Symmetry along the cuts. If £ belongs to the rays passing through the square-root branch points, 
(|22b[) implies that (|2Tj) is a polynomial equation featuring real coefficients. Therefore, its roots are 
real or in complex conjugate pairs. Consider the semiline cuts 7j, j = l,...,q, defined by 

(23) 7 , = U, arg£ = arg£«, |£| > |^|} , j = 1, ...,q. 

Then, if £ G jj, j = 1, q, two of the q roots lie on the segment (0, 1). If, in addition, £ sa 00, then 
one of these two roots belongs to the small circle around the origin and the other one belongs to the 
small circle around 1 (see Figure [TJ where the 00- configuration of the q roots is shown for £ s=s 00 and 
£671). 

3.3. Roots dynamics and topological properties. In order to understand the topological properties of 
the Riemann surface T, we first cut the ^-complex plane along the rays 7,-, j — 1, ...,q, see (f23|) . from the 
square root branch points £^ to the branch point at 00, and we introduce the closed contour 7 whose interior 
is the cut f-plane thereby obtained (see Figure [2]) . Let Fj be the Riemann sheet of Y associated to the root 
Wj. Our goal is to construct the q images Xj, j — 1, ...,<?, of the cut £-plane, corresponding to the q roots 
wj(£) of ([2"T]) . j = 1, q, and to study their connections. 

Notation: Here and hereafter, through Sections 13.21 and 15.21 all the indices in the designation of the square- 
root branch points ^ , the roots Wj , the images Xj and the sheets Tj are defined m5d (q) via the congruence 
defined in (J2J, namely = ^ MOD (<?)) or T 3 = 7} m6d 

Our strategy will be to start from the 00- configuration, namely the root configuration entailed by setting 
£ ~ 00, and to map out the structure of the Riemann surface as £ is moved in along rays and it is moved 
around along circles. 

Starting with £ 6 jj, it is natural to call Wj the root lying on the segment (0, 1) and belonging to the 
small circle around 1. Then Wj+i, w,j +q ^ p -i are the other roots of this circle, enumerated sequentially 
in counterclockwise order. Analogously, we denote by w.j+ q - p the root lying on the segment (0, 1) and 
belonging to the small circle around 0; and by Wj +g _ p+ i, Wj +q ^i the other roots of this circle, enumerated 
sequentially in counterclockwise order (see Figure [l}. 
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Using the large £ asymptotics, (|15al) - (|15cl) and (|16al) - (|16cl) . and the symmetries of T described above, we 
infer the following basic motions. 

(i) 

• As £ moves along the cut jj, from oo to £jr , the two roots uij and Wj+ q - p , lying on the segment 

(0, 1), move along it, from the small circles around and 1 to the collision point fi. In addition, being 

( i) 

£k a branch point of square- root type, a 2ir rotation of £ around it corresponds to a ir rotation of 
these two roots around \i. All this implies that, if £ travels along the contour surrounding the cut 
7j (for instance, moving from the point A to the point C along the path shown in Figure [2]), then 
the two roots Wj and Wj+ q - p involved in the collision exchange their position (see Figure [3]). The 
remaining roots are essentially unaffected by this motion, moving back and forth on lines and going 
back to their starting positions. We have thereby established the first basic motion: a motion of £ 
around the branch cut "fj yields an exchange of the two roots Wj and Wj+ q - p . 

• If, starting from the cut jj, £ performs a 2ir/q counterclockwise rotation, moving from the cut jj to 
the cut 7j+i, then the (q — p) roots surrounding 1 undergo a clockwise rotation around 1, while the p 
roots surrounding undergo a clockwise rotation around 0. When £, starting from 7^, reaches Jj+i, 
two new roots belonging to the two small circles get aligned on the segment (0, 1); they are just Wj+i 
and Wj+i+ q -p (see Figure [4]). Therefore the two sets of roots undergo cyclic permutations, which is 
the second basic motion: a rotation of £ from jj to Jj+i produces a cyclic permutation of the two 
sets of roots in the 00 -configuration {wi, w q - p } and {w q - p +i, w q }. 

Repeating q times the above two motions with respect to the other cuts, in sequential order, the point £ draws 
the whole closed contour 7 and correspondingly, due to the above-mentioned basic motions and symmetries, 
each root wj draws the closed contour in Figure [5] around the cut [0, 1] of the w-plane. The following details 
can be given (compare Figure [5] and Figure [SJ . 



"10 


"w 9 . ' W <> 

W 7 

w 5 


( 

w ll 


^ " ° ■ e r- 

W 4 

• 

W 

w 12 w 3 



Figure 3. Exchange of a pair of roots. As £ travels on the contour surrounding the cut 
7j (for example moving from the point A to the point C along the path shown in Figure 
HJ, Wj and Wj+ q -p interchange their positions, while the other (q — p) roots are essentially 
unaffected by this motion. 



(1) As £ moves on 7 through the points A, B and C, around the cut 77, Wj moves along the cut [0, 1] 
through the homologous points, exchanging its position with the root Wj+ q - p . 

(2) As £ moves on 7 from C to D, the relevant part of the motion of Wj consists in a clockwise rotation 
around 0, from C to D. 

(3) As £ moves on 7 through the points D, E and F, around the cut ^j- q + P , Wj moves along the cut 
[0, 1] through the homologous points, exchanging its position with the root Wj^ q+p . 

(4) As £ completes the contour 7, moving from F to the starting point A, also Wj completes its closed 
contour around [0, 1], and the relevant part of this motion consists in a clockwise rotation around 0, 
from F to A. 

From the above considerations we finally infer the following 
Topological properties of T for /i = p/q and < p < q. The Riemann surface T defined in S21\) is a 
q-sheeted covering of the £-plane of genus 0. In the finite part of T there are q square-root branch points: 
(£p , fi) € r, j = 1, ...,q (see \2S\) ). If Fj is the sheet associated to the root Wj, then the j-th branch point 
connects the sheets Tj and Fj+ q - p . Each sheet Fj contains just the two square-root branch points , fi) 
and (£p q+p \[i). The compactification ofT is achieved at £ = oo, where the branch point (00, 1), of order 
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Figure 4. Cyclic permutation of the two groups of roots. If £ £ jj, the roots Wj and Wj+ q ~ p 
are aligned on the segment (0, 1); after a 2-7r/g counterclockwise rotation of £, from the cut 
7j to the cut 7j+i, the two groups of roots {wj, .., Wj+ g _ p _i} and {tyj+ g _ p _i, .., 
undergo a clockwise rotation. When £ reaches 7j+i, then the roots Wj+i and w^+i+q-p get 
aligned on the segment (0, 1). 




Figure 5. As £ travels on the closed contour 7 in Figure [H each root travels on a closed 
contour around the cut [0, 1] in the iu-plane. Therefore the image of the cut £-plane is the 
w-plane cut along the segment [0, 1]). 



(q — p — 1), connects the first (q — p) sheets, and where the branch point (00, 0), of order p — I, connects the 
remaining p sheets. As £ turns counterclockwise around (00, 1), the connected sheets are visited in the order: 
J-j, J~j+ P , J-j+ip,...; as £ turns counterclockwise around (00, 0) instead, the connected sheets are visited in the 
order: J-j, J- j+ q - p , J~ j+2(q-p)>--- 

3.4. The physical root. Once a labeling system has been introduced for the square-root branch points 
£b s and for the roots «y's, in order to solve our original problem (JTJ we need to indicate which of the q 
roots of ([21]) corresponds to the physical root w(t) = w(£) appearing in (|6a |) - (|6bj) and satisfying (jSJ. 

As implied by (|12ap and (|12b|) . at t = we have that £ = £0 and one of the q roots of (j2"Tj) is wq. Therefore, 
it is sufficient to label wq in order to label the physical root. Observing (|12a[) and (|9b[) . we immediately note 
that £0 is a 9" valued function of the initial conditions zi(0), z 2 (0) and 23(C)); on the contrary from (|12b[) we 
infer that u> is a single- valued function of the initial conditions zi(0), z 2 (0) an d 23(0). Via (|5c| and (|10b[) . 
we get that £ an d £, the center of the circle 5 on which the £ variable moves, have the same degeneracy. 
This degeneracy can be removed by arbitrarily fixing one of the q determinations of R in (|9b[) , placing the 
circle S on the £-plane and then labeling counterclockwise the square-root branch points lying on the circle 
B defined by (|22b|) so that (f^ is the first branch point encountered by moving counterclockwise on the arc 
of the circle B obtained by intersecting the two circles and contained inside the circle 2 (see Figure [5]). When 
the two circles S and B do not intersect, then the first branch point £^ on B can be chosen arbitrarily. 

To proceed in labeling the physical root, we rephrase (|2ip in the following form 

(24) [(x — l) 2 + y 2 ] 3 ~* L (x 2 + y 2 )^ e i ' 9 arg ^ +( - q ~ p ^ arg i x ~ 1 + i v)+P ar g {x-Hy)] _ ^ ^ 

where x — 5f (w) and y = $s(w) are respectively the real and imaginary part of w and arg (£) is the argument 
of the complex variable £. Equation (p4|) implies that, for fixed £ and for rational values of /i, all the roots 
of (|2ip lie on the algebraic curve on the w-plane given in implicit form by 



(25) [(.-i) 2 + y 2 r p (, 2 + y 2 ) p Hr 29 
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Figure 6. The variable £ moves counterclockwise on the circle S starting from £q- The first 
square-root branch point on the circle B is the first one encountered by moving counter- 
clockwise on the arc of B, obtained by intersecting the two circles and contained inside the 
circle S. In this example p = 5 and q = 12. 

As implied by the discussion in Section [3721 for |£| < r& the curve (|25|) consists of a single closed branch that 
becomes a big circle of radius o(|£| _1 ) around the origin of the w-plane as |£| — > 0; for |£| > r& the curve (|25|) 
consists of two closed branches that become two small circles, one of radius 0(|£| -p ) around the origin and 
the other of radius 0(|£| _ «^) around 1, as |£| — > oo (see Figure [7]). 




(c) |£| = r b (d) |f | > r b 



Figure 7. The locus of the roots (j2l))) on the w-plane for different values of |£| when p = 5 
and q = 12. 

If the independent variable £ moves counterclockwise on a circle around the origin, namely if only the phase 
of £ varies, then the roots of (|2"TT) move clockwise on the curve ([2"51) . undergoing cyclic permutations: if 
(|25p is constituted by two unconnected branches, then the two systems of roots undergo two separate cyclic 
permutations. 

On the other hand, (|22bj) implies that, if £ moves along one of the rays that start at the origin of the 
£-plane and intersect the branch points, i.e. if arg(£) = arg(£^) for some j, then £ q e R. Therefore, via 
(pM|) . one immediately gets that in this case all the roots of (|2"TT) lie on the w-plane on the locus defined in 
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implicit form by 



(26) 



with 



(q — p) arg (x — 1 + i y) + p arg {x + i y) = nir 
1=1,3, 5, 2 q — 1 if (q — p) is odd 



n = 0, 2, 4, 2 <? — 2 if (g — p) is even . 

The branches of the curve (f26| define a partition of the w-plane (see Figure [8]). In other words, the partition 
of the £-plane into the q sectors obtained by tracing the rays that start at the origin and cross the square-root 
branch points (see Figure 9(a) I induces a partition of the ui-plane into the q sectors individuated by (1261) : if 



the variable £ moves entirely in a single angular sector on the £-plane, without crossing the lines that define 
the sectors, then each root of (|21[) is contained in a single region of the w-plane, without crossing any of the 
branches of the curve (l26t that separate the regions. 




Figure 8. The curve (|26[) divides the w-plane in q sectors. For a fixed value of £, each root 
of (f2~Tj) lies inside a sector. When £ moves counterclockwise around the cut jj, the roots Wj 
and Wq-p+j (that are respectively in the lower and upper exchange region) exchange their 
positions in the corresponding sectors. In this example p = 5 and q = 12. 

For q > 3 there is always a branch of the curve (f2"6"]l that intersects the point /j, on the w-plane. As 
shown in the previous section, when f moves around the branch cut then the two roots Wj and Wj+ q - p 
exchange their positions moving around the point /i on the w-plane (see Figure [8]). Therefore, if the root Wj 
is contained in the lower half-plane of the w-plane in the region that is bounded by the real segment (/x, 1) 
and by the lower branch of the curve (|26|) that crosses the point \i (the lower exchange region), then the root 



2-p is contained in the upper half-plane of the w-plane in the region that is bounded by the real segment 



(0,/i) and by the upper branch of the curve (|2l)|) that crosses the point \i (the upper exchange region). 

We assign to each angular sector on the £-plane the same label of the branch point that is next to the sector 

(counterclockwise), see Figure 9(a) If the initial datum £ is contained in the sector j on the £-plane, then 



correspondingly the root Wj is in the lower exchange region on the w-plane. Once a label has been assigned 
to the lower exchange region, we can label all the remaining regions (and correspondingly the roots in them) 
by counting them counterclockwise mod (q), see Figure [9(b)] When all the regions on the w-plane have 
been labeled, then the label of the physical root is the label of the region where wq falls. 
The prescription to label the physical root can be summarized as follows (see Figure [9]): 

(1) Trace the circle B in the £-plane, on which the branch points lie (see ([M)) ). 

(2) Choose arbitrarily the same determination for £o and R (see (|12ap ). 

(3) Trace the circle S (see (|10b[) ) and select, as explained in the present section, the first branch point 

■ Starting from there, label the remaining branch points counterclockwise. Define the angular 
sectors on the £-plane. 

(4) Compute £o and put it on the £-plane. Denote by j the sector where £o falls. 
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(a) Sector labeling of the £-plane. 



(b) Corresponding sector labeling of the w- 
plane. 



Figure 9. The enumeration scheme for the roots of (1211) . In this example p = 5 and q = 12, 
£o is in the fourth sector on the £-plane, while the physical root wq lies in the eleventh sector 
on the w-plane. 



(5) Trace the curve (|26p on the w-plane. Enumerate the regions defined by the branches of the curve by 
assigning m5d (q) the label j to the lower exchange region and moving counterclockwise. 

(6) From the initial data, compute Wq and put it on the w-plane. The label of the physical root is the 
label of the region where wq falls. 

This mechanism to identify the physical zero cannot be applied when fi is irrational as the separatrices of 
the sectors become infinitely close. Therefore this case requires a separate treatment. 



4. Connection with graph theory 

In this section we apply graph theory in order to describe the monodromy group associated to (|2ip , namely 
the subgroup of the symmetrical group S q corresponding to all possible exchanges and cyclic permutations 
occurring on the system of q roots when the independent variable £ moves along a closed path on the £-plane. 
Since eventually £ shall travel along the circle 3 described by (|10b[) and the square-root branch points lie 
on the circle B, see (|22c[) . we only need to consider the inclusions of (one or more) consecutive square-root 
branch points on B. 

We associate a labeled planar graph V to each mutual configuration of the two circles 3 and B. In order 
to construct this graph, let us trace q nodes - corresponding to the sheets of the Riemann surface T - on the 
vertices of a regular polygon; let us arbitrarily give to one of the nodes the label 1; if Vj is the label of the 
j-th node of the graph, counting counterclockwise from the node 1, then the labels are given by the following 
rule: v± = 1 and fj+i = (Vj +q—p) mod (g), j = 1, ...,(?— 1. When 3 includes no square-root branch points, 



then the planar graph has no edges (see Figure 10(a) I. If the circle 3 includes the square-root branch point 
$ , j = 1, ...,q, correspondingly the planar graph gets an edge between the nodes j and (j + q — p) mod (q) , 
as illustrated in Figure [TU1 With this prescription, V depends on 3, p and q, namely V = V(S,p, q). 

From the general theory we know that the inclusion of exactly q square-root branch points inside the circle 
3 is equivalent to the inclusion of the branch point at infinity: in this case the period of the zeros of (|2ip are 
q—p or p, see (|15a[) - (|15c[) respectively (|16a[) - (|16c|) . with the first q — p roots having period q — p and the last 
p roots having period p. From now on, we will suppose to have at most q — 1 consecutive square-root branch 
points, starting from ^ , inside the circle 3. 

A path of length n on the graph V is a sequence of n adjacent distinct nodes {vj 1 ,Vj 2 , ...,Vj n }, namely 
such that {vj 11 Vj 2 ), (vj 2 ,Vj 3 ), (vj n _ 1 ,Vj n ) are edges of the graph. For fixed 5, p and q, if in V(S,p,q) the 
node Vj belongs to a path of length b, then the period of the corresponding root of the algebraic equation 
?Tj) , w Vj (£), is exactly 6, while £ moves along the circle 3. If the node Vj is unconnected (i.e. it is touched 
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Figure 10. Planar graphs associated to different inclusions of consecutive square-root 
branch points in the circle S, for p = 5 and q = 12. In each caption, the included branch 
points are indicated. 



by no edges), then the period of w Vj (£) is 1. For fixed p and q and for a certain mutual configuration of the 
two circles EE and B, in order to describe the behavior of the roots of (|2"Tj) . we need to measure the lengths 
of the paths in the corresponding graph V(£,p, q). 

For a graph ^(EE,^, q) with q nodes, the set of all path lengths corresponds to a decomposition in cycles of 
elements of the symmetrical group S q . Once the correspondence between the graphs V"s and the elements of 
S q has been established, we can use the well-known correspondence between Ferrer diagrams and permuta- 
tions. We recall that a Ferrer diagram is an ordered disposition of q blank boxes in rows and columns, such 
that the number of boxes in each column equals the lengths of the cycles of the corresponding permutation 
and the total number of columns is the number of cycles in which the permutation is decomposed. For 
instance, the Ferrer diagram 



14 



D. GOMEZ-ULLATE, P. M. SANTINI, M. SOMMACAL, AND F. CALOGERO 



corresponds to the planar graph presented in Figure 10(g) In a Ferrer diagram, the columns are sorted in 



decreasing order of length, from left to right. Moreover, note that each column in a Ferrer diagram is oriented 
up-down, the first row being at the top, the last row being at the bottom. Then the problem of knowing the 
period of a root of equation (|2~Tj) . while £ moves along the circle 3, is reduced to measuring the lengths of 
the columns of an appropriate Ferrer diagram. 

4.1. The bumping rule. In this subsection we show how to recursively build the Ferrer diagrams associated 
to the planar graphs Vs. For fixed p and q, suppose that the circle H includes b square-root branch points, so 
that V(3, p, q) is the corresponding labeled planar graph. In the following, we include only consecutive square- 
root branch points of B in S, always starting from and moving counterclockwise: {£b , > •■•> £b }• 
Under this hypothesis, V(H,jj, q) depends on p and q and on 3 only via the number b of included square- 
root branch points. Assuming that p and q have been fixed once for all, we set V(E,p, q) = V(b). Let F(b) 
indicate the Ferrer diagram corresponding to the planar graph V(b). In what follows, we construct a recursive 
bumping rule for Ferrer diagrams, namely a rule that permits to build F(b + 1) from F(b). 
First of all, we need two integer sequences: 

(27a) q k = p k _ x , q = q , 

(27b) pk = qk-l MOD (pfc-l) , po — q - p , 

with < k < k, where k is an integer number such that 
(27c) q k = 1 ; 

note that k always exists due to the decreasing nature of the sequence {qk}- We need also an auxiliary 
recursive sequence, written as a combination of the previous two, 

(28a) b k = b k -i + qk-l - Pk-i , b = . 

We remark here that the 6-sequence (|28a|) divides the discrete segment [1, q— 1] into k parts of length qk—pk- 
Moreover, from the first of the (|27a|) . via (|28a[) . one gets 

b = o , b k = b k -i + qk-i - qk = b k -2 + qk-2 - qk-i + qk-i - qk 
= b k -2 + qk-2 - qk = — = h-h + qk-h - qk 

(28b) = b + q - q k = q - q k ■ 

Suppose that 3 includes b square-root branch points. If we modify 5 so that it includes also the next 
adjacent branch point, correspondingly we modify V(b) into V(b + 1), adding to V(b) a new edge between 
the nodes (b + 1) and (b + l + q — p) mod (q). 

If < 6 < 6i — 1 = q — p— 1, then each new edge added to the graph V(b) connects a single unconnected 
node to a path composed by a certain number of edges (at least, to another single unconnected node). When 
b = bi — 1, it is impossible to trace a new edge on V(b) connecting two unconnected nodes or a single node 
to a path: the insertion of a new edge in a graph at this point causes the connection of two paths, resulting 
in a new longer path in the graph V(b + 1). 

To understand what happens when b\ < b < 62 — 1 5 let us imagine to build an auxiliary graph in the 
following way: associate a weighted planar graph V(b) to the graph V(b), so that, for each path on V(b), 
you have a weighted node in V(b), with weights on V(b) equal to the lengths of the paths on V(b) (we recall 
that a single node is equivalent to a path of length 1). Now the discussion of the case b\ < b < 62 — 1 is 
analogous to the discussion of the previous case < b < bi — 1, using V{b) instead of V(b), with the only 
difference that the length of a path in V(b) is not just the number of nodes touched by the path, but the 
sum of all the weights of the nodes touched by the path. If b\ < b < 62 — lj then each new edge added to the 
auxiliary graph V(b) connects a single unconnected weighted node to a weighted path composed by a certain 
number of edges (at least, to another single unconnected weighted node). But when b — 62 — 1, again it is 
impossible to trace a new edge connecting two unconnected weighted nodes or a single weighted node to a 
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weighted path on V(b): the insertion of a new edge in V(b) at this point will cause the connection of two 
weighted paths, producing a new longer path in the auxiliary graph. 

In order to progressively include more branch points inside S, we must build each time a new auxiliary 
graph, in the afore-described way. We iterate this operation, building a new auxiliary weighted graph every 
time we arrive to include exactly bk branch points inside the circle 5, until we reach b = q — 1. In this 
language, qk is the number of unconnected weighted nodes in each auxiliary graph when bk < b < frfe+i, while 
Pk is the number of weighted nodes in each auxiliary graph when exactly bk+i branch points are included in 

The use of Ferrer diagrams strongly simplifies this picture. When b — bk, the lengths of the columns of the 
Ferrer diagram F(b) represent the values of the weights of the nodes of the corresponding auxiliary weighted 
graph. If bk < b < bk+i, then the lengths of the columns of the corresponding Ferrer diagrams F(b) represent 
the lengths of the weighted paths in the corresponding auxiliary weighted graphs. When b = bk, the number 
of columns in F(b) is qk- Moreover, note that, if one starts to count columns from the left in F(b), when 
bh < b < bh+i, then the columns lying at the positions ph + 1, Ph + 2,..., have all the same length, since they 
coincide with the shortest paths on the corresponding planar graphs. 

If the circle S includes b consecutive square-root branch points, when bj < b < fej+i, and if we modify 
it in order to include 6+1 consecutive square-root branch points, then we can build F(b + 1) just moving 
the whole last column on the right of the Ferrer diagram F(b) under the first available column from the left, 
never occupying positions beyond the Pj-th. column. To pass from bj to see (|28a|) . we have to move 

Qj ~ Pi columns. All the above considerations lead to the following prescription. 

The bumping rule. The Ferrer diagram F(0), corresponding to the inclusion of zero branch points, is 
composed by a single row of q columns of length 1. The Ferrer diagram F(b), with b < q—1, can be obtained 
from the Ferrer diagram F(b — 1) in the following way: 

• compute h such that bh <b < bh+i', 

• counting the columns of F(b — 1) from the left, move its whole column located at the position ph + 1 
under the column located at the position (b — bh) MOD (ph)- 

We remark that, in principle, one could ask to move any column lying at the position ph + 1, Ph + 2, since 
they all have the same length. The reason for prescribing a movement of the column lying at the position 
Ph + 1 will be clear after the introduction in Subsection 15 . 1 1 of a numeration for the blank boxes of the Ferrer 
diagrams. 

For each fc between and fc, ()27cj) . the Ferrer diagram F(bk) features qk columns of two possible lengths 
only, pk columns of length 7^ and qk — Pk columns of length T^ (see the diagrams in Figure [TTJ neglecting 
the numeration, that will be explained below). A k-level Ferrer diagram is a Ferrer diagram with exactly qk 
columns. 

The bumping rule implies that we "pass" from a fc- level Ferrer diagram to a (A; + l)-level Ferrer diagram 
once we have moved all the shortest columns on the right under the tallest columns on the left. 

The column lengths, T [ k 1] and T^\ in a fc-level Ferrer diagram can be given by the following recursive 



rule: 



(29) 



,7,(1) . ,7,(2), x ^(1) 1 

T k Pk + T k \1k ~ Pk) = qo, To = 1 



_ 71(1) , T (2) 
± k+l — ± k "r" ^k 



I ) . T (2) = 1 



where [x\ is the floor of the number x (namely, the largest integer less than or equal to x). Indeed, the first 
of the (|29p is a sort of conservation rule for the number of boxes in a fc-level Ferrer diagram. The second of 
the (|29p comes directly from the bumping rule: at each level, the length of the first column is the sum of the 
length of the column at the previous level plus the length of the last column multiplied by the number of the 
necessary column movements. Although relations (|29p are enough to determine the two quantities T^> and 
T^ 2 \ it is more convenient to put them in a different form: 



(30) 



_ , 7,(2) 

A k+1 ~ ± k "r J fc 



71(2) _ T (l) , T (2) 
± k+l ~ ± k "r ^k 



T y ' — 1 



T 



(i) 

(2) 
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Figure 11. Numbered Ferrer diagrams associated to different inclusions of consecutive 
square-root branch points in the circle S, for p = 5 and q = 12. In each caption, the 
number of included branch points b is indicated. The numeration is explained in the text. 

In the next subsection, we present a useful link between these two-steps-recursive relations, and the continued- 
fraction expansion of the number j^- . 

4.2. Connection with continued fractions. Let us recall some basic notions about simple continued 
fractions. Let x be a non-negative real number. We associate to the real number x an integer sequence {a^} 
such that: 

(31) x = ciq 



01 



(12 + 7 



»8 + - 

We say that (do, oi, a?,, ■•■) is the simple continued fraction of x with positive integers Vfc > (as usual, for 
the sake of simplicity, we avoid repeating the adjective simple when referring to simple continued fractions). 
The elements of the continued fraction expansion are called partial quotients. 

The integer sequence {a^} is finite if and only if the corresponding number a; is a rational number. The 
number of elements in the sequence is called the length of the continued fraction. 

The partial quotients can be recursively calculated introducing the auxiliary sequence of the remainders 



(32) 



r k — r , n j r o — 



i"fc-i-L>"fc-iJ 
ak = [r k \ 



We trust no reader will be confused by the similarity of this standard notation for the remainders rk with 
the radius r;,, see (|22cp . The rational number c„ = c n (x) obtained truncating the continued fraction of x at 
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the n-th term is called the n-th convergent of the continued fraction: 



(33) 



Pn 



c„ = — — [ao,ai, ...,«„) = a H 



The numerator P n and the denominator Q n of the n-th convergent c n satisfy the following second-order 
recurrence relations: 

^g^-j Pn = O-nPji-l + Pn-2, P-2 = 0, P-l = l] 

Qn = dnQn-1 + Qn-2, Q -2 = 1, Q-l = 0. 

Now we are ready to rephrase the quantities introduced in the previous subsection in the language of 
continued fractions. First of all, let us eliminate from the g-sequence the dependence on the p-sequence in 
formula (|27a|) : 

(35) q n = q n - 2 MOD (g n _i), q = q, qi=q~P- 

We can also write this formula as follows: 



(36) 

where we used the fact that 
(37) 



'In 



q n -2 



<?n-l Qn-1 

x-[x MOD (y)] 
V 



q n -2 

Qn-l 



qo = q, qi = q - v , 



| J for all x, y G N+. Now, let us set 



Qk+i 



with < k < k 



Using (13^1) . we recover formula (|32p for the quantity r^, 

1 



<Zo q 
witn rp = — = 



1 



(38) r k = - 

rh-i - \Tk-i\ 9i q-p 1-A* 

The q-sequence, namely the number of columns in the Ferrer diagrams, can be obtained by inverting formula 



(39) 



qk 



r k -i 



i 



T* 



qo = q, r 



i-m' 



< k < k 



where the r/~ are the terms of the sequence of the remainders (see (|38[0 of the continued fraction expansion 
of the number j^— , with partial quotients 



(40) 



a-k = \r k \ 



qk 
Pk 



and length k. The p-sequence can be calculated by using the relation: 

(41) p k = qk+i , < k < k . 

We can moreover provide a convenient reformulation yielding the lengths of the columns in the Ferrer 
fandTf 

(42) 

and the insertion of this last formula in the second of the ([30]) yields, via (|40|) . a 3-terms linear recursion for 

T-(2). 
± k ■ 



diagrams (13"0) . and Indeed, by subtracting the two recursion relations (|3"u) . one gets 

T (l) _ T (2) T (2) 
J fe+1 — ± k+l ' J fe ' 



(43) 



T (2) _ „ T (2) , T (2) 



A comparison of this last recursion relation with p4p . and an analogous comparison of the starting conditions, 
yield the relation 



(44) 



T { k 2) = Pk-i , 



where Pk is the numerator of the fc-th convergent of the continued fraction expansion of jrr;- Then, (J42 
and (UU imply 



(45) 



T, W - Pfc-l + Pk-2 
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Now all the quantities involved in the description of the Ferrer diagrams (i.e. the numbers of columns and 

boxes in each column) are expressed in terms of the continued fraction expansion of the (rational) number 
1 

1-/*' _ 

We conclude this section by noting a more convenient way of writing the fundamental 6-sequence (|28a[) . 
Indeed, by combining ([37]). (l38l) and (|40|) , one gets: 

(46a) q k = qk-2 - a k - 2 qk-i , qo = Q , Qi = q — p, 

and by comparing this last relation with the recursion relations (|34[) via the following ansatz: 

(46b) q k = {-if (aP k -2 + PQk-2) , k > , 

one immediately gets 

(46c) a=(p-q), P = q; 

the insertion of (I46b[) and (|46cl) in (I28b[) entails 

(47) b k =q-q k =q- {-lf[(p-q)P k -2+qQk-2} , k>0. 

5. Determination of the period 

5.1. The period formula for /i € Q. Once the structure of the Riemann surface has been achieved, our 
next step is to provide an explicit formula to predict the period of each of the q roots of the algebraic equation 

£3- 

The variable £ moves along the circle 5 which includes b consecutive square-root branch points on the 
£-plane, starting from Now we introduce a numbering of the boxes of the Ferrer diagrams, so that, if 
one fixes a single column, all the roots labeled with the numbers appearing in the chosen column have the 



same period, equal to the length of the column. For instance, look at Figure 1 1 (i) when 6 = 8. We interpret 
this picture, inferring that the roots w±, W3, wq, w% and wiq have period 5; the roots W2, wy and wg have 
period 3; the roots W4, W5, u>n and w±2 have period 2. This is exactly what you obtain by comparing with 



Figure 10(i) 



For fixed p and q, if we enumerate the boxes of F(0), starting from the left, with the first consecutive q nat- 
ural numbers (l,2,...,g) and if we apply recursively the bumping rule, then we obtain the correct enumeration 
for each F(b) with b > 0. 

Let us introduce the number lh(s), defined recursively as follows: 

(48) lh{ s ) — lh — l( s ) MOD (ph — l) , l (s) = 8, 

for all s = 1, 2, q and for all h = 1, 2, k. This number lh(s) gives the position of the column where the 
box with the number s lies in the /i-level Ferrer diagram. Once we know this information for a /i-level Ferrer 
diagram, we can extend the same information for a generic Ferrer diagram. 

Let us suppose that there are b branch points included in the circle S, with bh < b < b} l+1 . We are 
interested in the period of w s . Then two cases are possible. 

( 2) 

Case 1: If lh(s) > b — bh + pn then ' is the length of the column where the box with number s 
lies; this column was not moved from the position it occupied when exactly b^ branch points were 
included in the circle S. 

Case 2: If lh(s) < b — bh + Ph then we must consider two subcases: if lh(s) < Ph, then the length 
of the column in F(ph), containing the box s and lying at the position lh(s), is T^; if lh(s) > Ph, 
then the length of the column in F(bh), containing the box s and lying at the position ift(s), is 
( 2} 

T^>. In both cases, in order to obtain F(b) starting from F(bh), the bumping rule predicts that, 
if the circle S includes other b — bh branch points, then some of the last (qh — Ph) columns in 
F(bh), with length , are moved under the column located at the position lh(s) MOD (pft). If 
lh(s) < Phi then lh(s) m5d (ph) = lh(s) and the column at the position lh{s) in F(b) contains the 
box s. If lh(s) > ph, then the column in F(bh), at the position lh(s) and containing the box s, 
is moved under the column in F(b) lying at the position lh(s) MOD (ph). Thus, in both cases, the 
column in F(b) at the position lh(s) mod (ph) contains the box s. On the other hand, the number 
of moved columns depends on whether lh(s) mod (ph) is smaller than (b — bh) MOD (ph) or not. If 

~ ~ (2} 

lh(s) MOD (ph) < (b — bh) mod (ph) then the number of columns of length T ?i , moved under the 
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column at the position l h (s) MOD (p h ), is j). If l h (s) mod (p h ) > (b — bh) mod (p h ) then 

one needs to move one more column of length so their total number is il b ^ b p h h ^ 1 J + !)• 
From the above considerations, we finally infer the following 

Theorem 1. Assume < ji < 1 and \i 6 Q. Let the roots of the algebraic equation H21\) be labeled as 
described in Subsection \3.3i Let T(s,b) be the period ofw s , the s-th root of the algebraic equation i21\) when 
£ moves along a closed path on the ^-plane, including b consecutive adjacent branch points starting from the 
branch point ■ Let h be the integer such that < bh < b < bh+i < q — 1 and lh(s) the symbol Let 
TP'' and T/ 2 ^ be the quantities described by the recursions H30\) . Then we have the following period formula, 
for all b < q: 



(49a) 



T(s,b) 



T, 



(i) 



b-b h 



T, 



.(1) 



b-b h -l 



T, 



(2) 



iflh(s) <b — bh+Ph and 
(b-b h ) MOD ( Ph ) < l h (s) mod (p h ) 



Ph 



1 ) if l h (s) < b - b h + p h and 



T-, 



,(2) 



If b = q we have 
(49b) 



(b-b h ) m5d (p h ) > l h (s) MOD (ph) 
iflh(s) > b-bh+Ph ■ 



{q-p ifl<s<q-p 
p if q — p + l < s < q . 

Via the connection with continued fractions, as shown in Subsection \4-2\ we can reformulate J^ffa| ) as follows: 

' P h - 2 + ([!^\+l) Ph-i *fk(s)<b 



(50) T(s, 



, „,. , _ J + q- (b h + b h+1 ) and 
(b - b h ) MOD (q - b h+1 ) < l h (s) MOD (q - b h+ i) 



Ph- 



h-2 



b-b h 



q-bh^ 



2 Ph-i 



if hi(s) < b + q - (b h + b h+1 ) and 
(b - b h ) MOD (q - b h+ i) > l h (s) MOD (q 



Ph-i if lh(s) >b + q-(b h + b h+1 ) , 

where Pk is the numerator of the k-th convergent of the continued fraction expansion of j—. 

For any inclusion of consecutive branch points, the final formula (|50p indicates that the roots of the algebraic 
equation (|21l) can feature only three possible periods (see Table 1). Note that the sum of the first and third 
expressions in the right hand side of this formula, (f5U|) . always gives the second expression. Moreover, if b takes 
the special values b = bh+n (q—bh+i) , < n < ah— I , n £ N, where ah is the h-th partial quotient appearing 
in the continued fraction expansion of ■=-=—, then the roots of (1211) can feature only two possible periods: indeed, 
in these particular cases, the first condition in (151)1) fails for all s since (6 — bh) MOD (q — bh+i) = q — bh+i 
and 1 < l h (s) MOD (q - b h+ i) <q - b h+1 . 

Summarizing, given the initial data Zj(0), i = 1,2, 3, and the coupling constants /, g (such that p <E Q) we 
can now predict the period of the function w(t) by the following prescription: 

(1) Draw the circles S and B defined by (|14bj) and (|10b|) . 

(2) Identify the label of the physical root s following the prescription given in Section I3T41 

(3) Calculate the number of branch points b included in the circle H. 

(4) Develop in continued fraction the ratio and build the sequence b^ defined by (|47|) . Identify the 
element h of the sequence such that bh < b < bh+i- Calculate lh(§) using l|4"8|) . 

(5) The period T(s, 6) of w(t) is given by Theorem 1. 

The understanding of the topology of the Riemann surface T for p 6 Q allows one to associate the evolution 
in time of the physical problem ([T]) to the symbolic dynamics given by the sequence of the labels of the visited 
sheets of T. In other words, one can describe the time-T map obtained by sampling the solution at every 
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Table 1. Values of T(s, 6) when p = 5 and q = 12. 



time interval T, as a sequence of natural numbers, i.e. the sequence of numbers labeling the visited sheets, 
after each round trip as £ travels on the Riemann surface. 



5.2. The period formula for p, ^ Q. We treat the case in which p is an irrational number as a limit of 
the case in which /i is a rational number. If /i = p/q gets close to being irrational, the number q - that 
corresponds to the number of roots of equation (|2ip as well as the number of branch points of in the 
complex £-plane - approaches infinity. So we must analyze all the previous formulae in the limit: 

(51) p, q — > oo with < fi = — < 1 , 

namely in the limit in which the integers p and q diverge so that their ratio tends to the irrational number \x. 

If p is an irrational number, it follows that the Riemann surface T associated to ([lip becomes an oo-sheeted 
covering of the complex £-plane. In the limit (1511) . the continued fraction expansion of yjj— has length k that 
tends to infinity: k — > oo. The first step is to normalize the g-sequence (|35|) . introducing the new sequence 

(52) p k = , po = 1 , pi = l-(i,, 

Tk-l 

where r k is the sequence of the (irrational) remainders of the continued fraction of j^- (see (1551) '). In the 
limit (I5ip , the p-sequence becomes a strictly decreasing sequence of irrational numbers such that p^ — > as 
k — > oo. 

The next step is to obtain the proper continuous variable which replaces the discrete index b. Let B be the 
circle on which the square-root branch points lie and 3 the circle on which the £ variable travels periodically, 
see (jlObj) . We define the variable v, which is if S does not contain any arc of B and 1 if 3 does contain the 
whole circle B. If B and 3 intersect, select the arc of the B circle which is contained inside the evolutionary 
circle 3 and denote by v the corresponding angle </>, normalized by 2-7T (see Figure I12p . Since 

(53) 

2n q 

this is a convenient continuous variable appropriate to replace b in the case p is an irrational number. 
We then normalize the fundamental 6-sequence (|T7|) . replacing it by the new sequence {vt} ~ {^j" 

(54a) v h = 1 - p k , ^o = 0, k e N ; 

in the limit (|51l) . we have that — > 1 as k — > oo. It is also convenient to reformulate v k using the language 
of continued fractions. From (|47p . in the limit (f5Tj) we have: 

(54b) ^ = l-(-l) fc [(/i-l)P fe _ 2 +Q fe _ 2 ] , fceN. 

We are finally ready to enunciate the period formula when p is an irrational number. 
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FIGURE 12. The definition of the angle 4> 



Theorem 2. Let < fj, < 1 and /i ^ Q. Let T(y) be the period of one of the roots of the algebraic equation 
ill]) when £ travels on the circle 3 in the complex £-plane, intersecting the square-root branch points circle 
B in such a way that < v < 1 (see H53\) ). Let h be the integer such that < < v < Vh+i < 1- Let be 



the numerator of the k-th convergent of the continued fraction expansion of ■ 
T(v) = 1. Otherwise, T(v) takes one of the following three values: 



(see (5$). Ifv = 0, then 



(55) 



P 



h-2 



P 



h-2 



1-Vh- 



1— "M 



1 P 



h-\i 



2 P 



h-i ■ 



P 



h-l ■ 



If v = 1, namely if the evolutionary circle 3 contains entirely the square-root branch points circle B, then 
the dynamics becomes again simple, because the evolutionary curve effectively surrounds only the irrational 
branch point at infinity. In this latter case, the time evolution of the generic root of ill]) is quasi-periodic, 
involving a (nonlinear) superposition of two periodic evolutions with noncongruent periods, 1 and — or 1 and 
— !— , as can be easily understood observing the exponents of the asymptotics U5a\) - [T5c]) and \16a]) - l[76c]) . 



Note that, as in the rational case, the sum of the first and third expressions in the right hand side of ([53]) 
always gives the second expression. We remark that in the irrational case there are no labels for the roots, so 
we cannot assign the periods to the roots in terms of their labels as we did in the rational case. Therefore, a 
specific prescription to identify the period of the physical root w(t) cannot be given. At best, we can assure 
that the period is one of the three possibilities specified by (|55|) . 

For v = 0, the period T(0) is equal to 1. Since 1 is the accumulation point of the sequence {vk}, if v 
is instead close to 1 (i.e. if the evolutionary circle 5 contains almost completely the branch points circle 
B), a small change in the time trajectory results in a drastic change of the observed periods. In fact, T(v) 
approaches infinity as v — > 1 and the values of T[v) depend on the partial quotients of the continued 
fraction expansion of yi^S such partial quotients are well-known to be chaotic and unpredictable in their 
sequence for a generic irrational number (for almost all the irrational numbers, except for the quadratic 
irrationals that feature periodic continued fraction expansions). 

In the next subsection we treat a remarkable example where it is possible to obtain explicitly the asymptotic 
behavior of T(v) as v ~ 1. 

5.3. A remarkable example. In this subsection we display the special example with the following (conve- 
niently chosen) quadratic irrational value of fi in the interval < fj, < 1: 

(56a) [i = 



such that 
(56b) 



1 



3 + VE i + y 

1 + V5 



= f = 



l-/i r 2 

where if is the so-called golden ratio, namely the positive solution of the second degree equation 



(56c) 



tp-l=0. 
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The golden ratio <p has the nice property that all the infinite partial quotients (|4"0")) appearing in its 
continued fraction expansion are equal to unity: 



(57) 



tp = 1 + 



1 



1 



a k = 1 Vfc > 0. 



So, in this particular case, the recursion relation (|38p becomes: 

1 



(58) r k = 

Vk-l - 1 

By combining this last relation with (|5"2")) we get: 



with ro = <p , k > 1 . 



(59) 

and from this, via (|54a|) . 
(60) 



Pk+i = Pk-i - Pk with po = 1 j Pi = — = <P ~ 1 , k>l 



v k +i = Vk-i - v k + 1 with v — 0, v x = \ , k>l. 



Solving this last relation with respect to k, we get 



(61) 



Vk = 1 - 



k > 0. 



Moreover - somewhat remarkably - ([57)) implies, via (|34|). the following relation for the numerators of the 
convergents: 

(62) P k = P k -! + P k -2 , P-2 = , P-i = 1 , k > , 

i.e. exactly the recurrence relation for Fibonacci's numbers. Using Binet formula, we have: 



(63) 



V k+2 {-<P) 



-(fc+2) 



, k > -2 



Through J61|) . we can explicitly invert the inequality v k < ^ < I'fe+i, finding, for a fixed value of v in the 
interval < v < 1, the integer number A: such that z/fc < v < Vk+i"- 



(64) 



jfe = k{v) = 



log(l-^) 
log((p) 



- 1 , if < v < 1 



, ifi/ = 0. 

Using (jMJ), we see that the argument of the floor function in (|55|) must satisfy the inequalities: 



0<z <<y9-l<l ; 



(65) 

J- - ' / fe( 1 y)+i 

so the floor function in (|55|) always vanishes and the root period T(v) has one of the following values: 

(66) T{y) = {P fe+ i, P fe , Pfe_i} with < i/ fc < v < v k +i < 1 , 

namely one of three consecutive Fibonacci's numbers. From (|62[) . via (l64l) . we see that, for < i/ < 1, 



(67) 



1 

7! 



l-v 



\-v 



< Pk{u) < 



V5 



l-v 



l-v 



From the above relations (1571) and (|65l) . and from the period formula (|55|) . we obtain the following lower 
and upper bounds for the period values in terms of v in the interval < v < 1: 



(68) 



v(2-v) <T , )< 7 + VE+(VE-3)v(2-v) 



V5(l-v) ~ 2V5(l-v) 
These inequalities entail that the integer T(v) diverges proportionally to (1 — v)^ 1 as v — > 1. 
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6. Summary and conclusions 

We have studied the trajectories of the system of three coupled ODEs (TIJ in the semi-symmetrical case 
((4]) when two coupling constants are equal, and the ratio \i defined in (0 belongs to the interval (0, 1). 

For rational values of fx, all orbits are periodic for arbitrary initial data, except for a set of null measure 
that corresponds to a collision in finite time. Moreover, the system is isochronous: the periodic orbits are 
neutrally stable and a small perturbation of the initial condition produces another periodic orbit of the same 
period. In this case, the Cauchy problem can be solved explicitly and the evolution in time of the physical 
trajectories can be put in correspondence with the symbolic dynamics (a sequence of natural numbers) of 
the labels of the visited sheets of the associated Ricmann surface. 

For irrational values of the ratio fi, one must distinguish two cases depending on the relative position of 
the evolutionary circle S, defined by (|10b|) . with respect to the branch-point circle B, defined by (|14b[) . If B 
is entirely contained in 5, then the solution is quasi-periodic, while if 3 contains only a part of B, then almost 
all orbits are periodic (again, except for a set of null measure that corresponds to a collision in finite time). 
The relative position of these two circles depends on the choice of initial data. In fact, initial data can always 
be found such that the period of the corresponding orbit is as high as desired (although always a multiple 
of the fundamental period) . This situation corresponds to the limiting case between the two cases described 
above. An explicit instance of this situation has been given for a particularly simple example related to the 
golden ratio and Fibonacci's numbers. 

The main result of this paper is to derive explicit formulae for the period as a function of the initial data 
and the coupling constants. The only way the period of an orbit can change as the initial data are varied 
is passing through a singularity (corresponding to a collision between the particles). The dependence of the 
period on the initial data is shown to depend on the continued fraction expansion of the parameter /i, which 
for irrational values shows a rich behavior. We would like to emphasize the novelty of the approach developed 
in this paper: it is not possible to achieve such an explicit description of the dynamics of the system using 
the qualitative theory of ODEs, which illustrates why the complete description of the Riemann surface is an 
important goal in the cases where it can be achieved. The description performed in this case involves the 
use of techniques that range from complex analysis and geometry to graph theory and continued fractions, 
which shows the novelty of the approach. 

One of the important lessons from the results in this paper in relation to the ideas of Kruskal, Bountis, 
Tabor and their collaborators mentioned in the Introduction is the following: the fact that the Riemann 
surface associated to the solution of the problem is infinitely sheeted is not enough for a system to produce 
chaotic behaviour, although it can be a source of unpredictable long term behaviour. This is illustrated by 
the behaviour of the semi-symmetric three-body problem with an irrational value of \x (i.e. generic values of 
the coupling constants). In this case, we have shown that the Riemann surface is infinitely-sheeted and it 
contains an infinite number of branch points whose position depends on the initial data and whose projection 
on the complex plane fills densely a circle B. Furthermore, we have proved that all solutions are either 
periodic, quasi-periodic or lead to a collision in finite time. Despite this fact, even in the periodic case there 
is a source of unpredictability associated to the fact that it is not possible to identify the root associated to 
the physical problem among the infinity of roots of dill) . Moreover, two trajectories originating from nearby 
initial data will separate at some future time as they will travel to different sheets on the Riemann surface, 
corresponding to the fact that the collision manifolds (locus of initial data in phase space that lead to a 
collision in finite time) fill densely an open set of phase space. On the other hand, in the /^-rational case, not 
only can the period be predicted by an explicit formula, but the complete motion can be described by the 
sequence of sheets in the Riemann surface that are visited during the motion. 

The description of the remaining case /i > 1 can be tackled with similar techniques to the ones developed 
in this paper. However, we have decided to postpone this analysis to a further publication, where the sensitive 
dependence on the initial conditions will be discussed. An interesting further step in this program would 
be to understand how sensitive dependence or classical indicators of chaos such as Lyapunov exponents are 
related to properties of the Riemann surface. 
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